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Motivated by a recent change in viewing the onset of the extra-galactic component in the cos¬ 
mic ray spectrum, we have fitted the observed data down to 10®'® GeV and have obtained the 
corresponding power emissivity. This transition energy is well below the threshold for resonant p'y 
absorption on the cosmic microwave background, and thus source evolution is an essential ingredient 
in the fitting procedure. Two-parameter fits in the spectral and redshift evolution indices show that 
a standard Fermi E~^ source spectrum is excluded at larger than 95% confidence level (CL). Armed 
with the primordial emissivity, we follow Waxman and Bahcall to derive the associated neutrino flux 
on the basis of optically thin sources. For pp interactions as the generating mechanism, the neutrino 
flux exceeds the AMANDA-BIO 90% CL upper limits. In the case of py dominance, the flux is 
consistent with AMANDA-BIO data. In the new scenario the source neutrino flux dominates over 
the cosmogenic flux at energies below 10® GeV. Should data from AMANDA-II prove consistent 
with the model, we show that IceCube can measure the characteristic power law of the neutrino 
spectrum, and thus provide a window on the source dynamics. 


I. INTRODUCTION 

A plethora of explanations have been proposed to ad¬ 
dress the production mechanism of ultra-high energy cos¬ 
mic rays [1]. In the absence of a single model which is 
consistent with all data, the origin of these particles re¬ 
mains a mystery. Clues to solve the mystery are not im¬ 
mediately forthcoming from the data, particularly since 
various experiments report mutually inconsistent results. 
In recent years, a somewhat confused picture regard¬ 
ing the energy spectrum and arrival direction distribu¬ 
tion has been emerging. Since 1998, the AGASA Col¬ 
laboration has consistently reported [2] a continuation 
of the spectrum beyond the expected Greisen-Zatsepin- 
Kuzmin (GZK) cutoff [3], which should arise at about 
GeV if cosmic ray sources are at cosmological 
distances. In contrast, the most recent results from 
HiRes [4] describe a spectrum which is consistent with the 
expected GZK feature. This situation exposes the chal¬ 
lenge posed by systematic errors in these types of mea¬ 
surements. Further confusing the issue, the AGASA Col¬ 
laboration reports observations of event clusters which 
have a chance probability smaller than 1% to arise from 
a random distribution [5], whereas the recent analysis 
reported by the HiRes Collaboration showed that their 
data are consistent with no clustering among the highest 
energy events [6]. 

Deepening the mystery, recent HiRes data have been 
interpreted as a change in cosmic ray composition, from 
heavy nuclei to protons, at ~ 10® GeV [7]. This is an 
order of magnitude lower in energy than the previous 
crossover deduced from the Fly’s Eye data [8]. The end¬ 
point of the galactic flux is expected to be dominated 
by iron, as the large charge Ze oi heavy nuclei reduces 
their Larmor radius (containment scales linearly with Z) 
and facilitates their acceleration to highest energy (again 


scaling linearly with Z). The dominance of nuclei in the 
high energy region of the Galactic flux carries the im¬ 
plication that any changeover to protons represents the 
onset of dominance by an extra-galactic component. The 
inference from this new HiRes data is therefore that the 
extra-galactic flux is beginning to dominate the Galac¬ 
tic flux already at ~ 10® GeV. Significantly, this is well 
below Eqzk ^ lO^®-"^ GeV [3], the threshold energy for 
resonant pycMB ^ > Nir energy-loss on the cosmic 

microwave background (CMB), and so samples sources 
even at large redshift. 

The dominance of extra-galactic protons at lower en¬ 
ergy can be consistent with recently corroborated struc¬ 
tures in the cosmic ray spectrum. A second knee, rec¬ 
ognized originally in AGASA data [9], is now confirmed 
by the HiRes-MIA Collaboration [10]. At 10®'® GeV, the 
energy spectrum steepens from E~^ to This steep¬ 

ening at the second knee can be explained [11] by energy 
losses of extra-galactic protons over cosmic distances, due 
to e+e" pair-production on the CMB. The theoretical 
threshold of the energy-loss feature occurs at 10®'® GeV, 
and therefore allows for proton dominance even below 
this energy. However, the HiRes data [7] seem to indi¬ 
cate a composition change coincident with the energy of 
the second knee (from about 50% protons just below to 
80% protons just above), and therefore argues for the 
beginning of extra-galactic proton dominance at the sec¬ 
ond knee. Another feature in the cosmic ray spectrum 
is the ankle at ~ 10^® GeV where the spectrum flattens 
from to This has been commonly identifled 

with the onset of the extra-galactic flux in the past. In 
the aftermath of the new HiRes data, the ankle can now 
be interpreted as the minimum in the e“'"e“ energy-loss 
feature. 

These changes in viewing the onset of the extra- 
galactic component have spurred a refitting of the cos- 
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mic ray data down to 10® ® GeV with appropriate 
propagation functions and extra-galactic injection spec¬ 
tra [11, 12]. The major result is that the injection spec¬ 
trum is signihcantly steeper than the standard E~'^ pre¬ 
dicted by Fermi engines. This result has consequences 
for neutrino observation: predictions for both the cos- 
mogenic fluxes (produced via interactions of super-GZK 
cosmic-rays on the GMB) and the direct neutrino lumi¬ 
nosity from optically thin sources can be signihcantly 
modihed. The implication for cosmogenic neutrinos has 
been discussed elsewhere [13, 14]. In this paper we an¬ 
alyze the impact on neutrino luminosities from optically 
thin sources which are associated with this change in view 
of the Galactic/extra-galactic crossover energy. 

The outline of the paper is as follows. We begin in 
Sec. II with an estimate of the power density of cos¬ 
mic rays assuming a “low” energy onset of dominance by 
an extra-galactic component. This is accomplished by a 
goodness-of-ht test of our scenario with the energy spec¬ 
trum as observed by the Akeno [9] -I- AGASA [2] and the 
Fly’s Eye [8, 15] + HiRes [4] experiments, respectively, 
in the energy range from the second knee upward, i.e., 
10® ® GeV < E < 10^^ GeV. In the htting procedure 
we use appropriate propagation functions [16] that take 
into account photo-meson and pair production on the 
GMB. Armed with the cosmic ray emissivities required 
to populate the observed spectrum with extra-galactic 
protons all the way down to 10® ® GeV, in Sec. Ill we 
derive predictions of neutrino fluxes associated with py 
or pp interactions in sources which are optically thin. To 
this end, we follow the procedure delineated by Wax- 
man and Bahcall (WB) [17], but instead of assuming a 
specific cosmic ray injection spectrum cx E~^ and red- 
shift source evolution cx (1 -I- z)®, we use values for the 
spectral and source evolution indices complying with the 
best fits to the spectra obtained in Sec. II. In Sec. IV 
we review the computation of cosmogenic neutrinos and 
show that in the “new” low crossover scenario, associ¬ 
ated neutrino spectra from optically thin sources domi¬ 
nates over the cosmogenic flux. In Sec. V we calculate 
event rates at IceCube [18] expected from source fluxes 
derived in Sec. Ill, and on this basis assess the potential 
of this detector to constrain the crossover energy between 
Galactic and extra-galactic dominance in the cosmic ray 
spectrum. Section VI contains a summary of our results 
and conclusions. 


II. EXTRA-GALACTIC COSMIC RAY POWER 

It is helpful to envision the cosmic ray engines as ma¬ 
chines where protons are accelerated and (possibly) per¬ 
manently confined by the magnetic fields of the accelera¬ 
tion region. The production of neutrons and charged pi- 
ons and subsequent decay produces both neutrinos and 
cosmic rays: the former via ^ (and the 

conjugate process), the latter via neutron diffusion from 
the region of the confined protons. If the neutrino- 
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FIG. 1: The goodness-of-fit test of the low crossover scenario, 
using the method outlined in Ref. [13]. In the left panel we 
show the 95% CL allowed regions in the n — 7 plane, obtained 
with Fly’s Eye -I- HiRes and Akeno -|- AGASA data. In the 
right panel we show the corresponding 68 % CL obtained from 
Fly’s Eye -I- HiRes data. In the fitting procedure we used 
data in the energy interval from E- = 10®'® GeV to E+ = 
10^^ GeV, taking Zmin = 0.012 and Zmax = 2. 


emitting source also produces high energy cosmic rays, 
then pion production must be the principal agent for the 
high energy cutoff on the proton spectrum. Conversely, 
since the protons must undergo sufficient acceleration, in¬ 
elastic pion production needs to be small below the cut¬ 
off energy; consequently, the plasma must be optically 
thin. Since the interaction time for protons is greatly 
increased over that of neutrons because of magnetic con¬ 
finement, the neutrons escape before interacting, and on 
decay give rise to the observed cosmic ray flux. The fore¬ 
going can be summarized as three conditions on the char¬ 
acteristic nucleon interaction time scale Tint; the neutron 
decay lifetime t„; the characteristic cycle time of con¬ 
finement Tcycie; and the total proton confinement time 
R:onf- (1) Rnt ^ R;yclei (2) > 'Tcyclej (3) Tint ^ T^onf ■ 

The first condition ensures that the protons attain suffi¬ 
cient energy. Conditions (1) and (2) allow the neutrons 
to escape the source before decaying. Condition (3) per¬ 
mits sufficient interaction to produce neutrons and neu¬ 
trinos. We take these three conditions together to dehne 
an optically thin source. A desirable property of this 
low-damping scenario is that a single source will produce 
cosmic rays with a smooth spectrum across a wide range 
of energy. 

Assigning extra-galactic dominance to energies begin¬ 
ning at ~ 10® GeV, rather than ~ 10®® GeV, increases 
the required energy production rate of extra-galactic 
cosmic rays. The power density in the energy range 
10®® GeV to 10®® GeV is found to be ecj^[10®®, 10®®] ^ 
5 X lO"®"® erg Mpc“® yr“® [19]. As emphasized in [19], this 
result is independent of source evolution: for the stated 
energy interval, cosmic rays from distant sources will un¬ 
dergo significant energy losses on the GMB, and thus only 
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nearby sources contribute to the observed spectrum. In 
what follows we obtain analogous power densities corre¬ 
sponding to the lower energy onset of extra-galactic dom¬ 
inance. As can be expected, these will have sensitivity to 
source evolution. 

We assume an isotropic distribution of neutron- 
emitting sources that can be described by a comoving 
luminosity distribution Cn{r,Ei) where Ei is the injec¬ 
tion energy and r the distance to Earth. The number of 
protons Np arriving at Earth with energy E per units of 
energy, area A, time t and solid angle reads, 


'A) — 




dE dA dt dil 
An 


dEi 


dr 


dPp\n{E-,Ei,r) 


dE 


( 1 ) 


where Pp|„ is the propagation function introduced in 
Ref. [16] which gives the expected number of protons 
above a threshold energy E if a neutron with energy Ei 
was emitted from a source at a distance r. The Monte- 
Carlo program that computes Pp\n uses SOPHIA [20] 
as an external package for simulation of the GZK in¬ 
teractions. The program exploits a continuous energy 
loss approximation to describe the e“''e“ pair production 
process. To estimate the differential flux of protons, we 
calculate the Pp\n function for infinitesimal steps (1 10 

kpc) as a function of the redshift z and multiply the corre¬ 
sponding infinitesimal probabilities starting at a distance 
r{z) down to Earth with z = 0. 

We take all sources to have identical injection spectra 
dNn/dEi cx 0 (ifi niax ~ Ei), and parametrize the 
redshift evolution of the source luminosity and the co¬ 
moving number density by a simple power-law. 


En — Pen [1 + ^(^)] Zmin) ©(-Zmax z) , (2) 

where the redshift z and the distance r are related by 
dz = (1 + z)H{z)dr. The Hubble expansion rate at 
a redshift z is related to the present one Hq through 
= Hq + z)^ + IIa] , where D,m and Ha are 

the matter and vacuum-energy densities in terms of the 
critical density. Here we take Hm = 0.3 and IIa = 0.7, 
in agreement with WMAP observations [21]. The results 
turn out to be rather insensitive to the precise values 
of the cosmological parameters within their uncertain¬ 
ties. The minimal and maximal redshift parameters Zmin 
and Zmax exclude the existence of nearby and early time 
sources. As default value, we take Zmin = 0 . 012 , corre¬ 
sponding to rmin = 50 Mpe, and comment on the effect 
of possible variations where appropriate. Note that the 
effects due to a change in Zmax can be largely compen¬ 
sated by a change in n. Therefore, we fix Zmax = 2 in 
the following and study the dependences on n only. For 
the maximum injection energy we take as a default value 
Ei^max = 10^^-® GeV. Most of our results are insensitive 
to this choice, as long as Ei^^ax is above ~ 10^^ ® GeV 
(see also Ref. [13]), as we will see in the following. 
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FIG. 2: Best fits to the ultra-high energy cosmic ray spec¬ 
trum in the energy interval [E-, i?+] as observed by Akeno 
-h AGASA and Fly’s Eye -b HiRes. We set E- = 10® ® GeV, 
E+ = lO" GeV, Zmin = 0.012, Zmax = 2, and Ei.max = 
10" ® GeV. 


The cosmic ray spectra as observed by Akeno [9] -|- 
AGASA [2] and Fly’s Eye [ 8 , 15] -I- HiRes [4] are fitted 
and confidence levels are assigned using a Poisson likeli¬ 
hood following the procedure detailed in Ref. [13]. The 
95% GL exclusion contours in the n — 7 plane for the two 
data samples are shown in Fig. 1. The parameters for the 
best fit, shown in Fig. 2, are given in Table 1. The dis¬ 
parity in the goodness-of-fit tests of the two data samples 
is largely originating in the presence of a spurious bump 
in the region of 10®'^ GeV of the Akeno -I- AGASA data 
(cf. Ref. [13]). This in turn stems from combining the 
data of the Akeno array and the full AGASA experiment. 
Interestingly, the new scenario with extra-galactic cosmic 
rays dominating the spectrum below the ankle, down to 
the second knee at = 10® ® GeV, is inconsistent at 
more than a 2a level with standard Fermi engine models 
that suggest an injection spectrum cx E~‘^. Note that for 
both data samples the best fit yields 7 = 2.54. Addi¬ 
tionally, 7 < 2.4 is disfavoured at the Icr level by Fly’s 
Eye [ 8 , 15] -b HiRes [4] data and at the 2a level by 
Akeno [9] -b AGASA [2] data. We have checked that 
this is robust against variations of Zmin < 0.012 and 
E^^n^ax > 10" ® GeV. 

It is conceivable that the extra-galactic proton flux 
largely exceeds the nearby data below 10®'® GeV. As can 
be seen in Fig. 2, this is not the case for the best-fit 
values. Moreover, if we make a more sophisticated anal¬ 
ysis with a Galactic component below 10® ® GeV, the 2a 
allowed regions will shrink. 
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TABLE I: Best fit parameters 


Experiment 

7 

n 

ecR[10® ® GeV, 10^^ ® GeV] 

AGASA 

2.54 

3.65 

2.5 X 10"'® erg Mpc“® yr“' 

HiRes 

2.54 

3.45 

1.3 X 10"'® erg Mpc“® yr“' 


At this stage, it is worthwhile to mention that we 
have verified the consistency of the simulations by fit¬ 
ting the data above 10^° GeV and comparing the cosmic 
ray power density, 

I’-Ei.ma.x 

^i,max] — PcR I dEi . (3) 

JEi^rnin dEi 

with the result obtained in [19]. Our best fits, us¬ 
ing Eq. (1) with 7 = 2 and n = 3, correspond to 
GeV, 10^^ GeV] = 4.8 x 10"^^ erg Mpc“^ yr“^ 
and GeV, 10^^ GeV] = 2.4 x erg Mpc“^ 

yr“^, for AGASA and HiRes, respectively. 


III. NEUTRINO COSMIC RAY 
CONNECTION 

For optically thin sources, the neutrino power den¬ 
sity scales linearly with the cosmic ray power density 
^cR The actual value of the neutrino flux depends 

on what fraction of the proton energy is converted to 
charged pions (which then decay to neutrinos). To quan¬ 
tify this, we follow WB and define as the ratio of 
charged pion energy to the emerging nucleon energy at 
the source. Depending on the relative ambient gas and 
photon densities, charged pion production proceeds ei¬ 
ther through inelastic pp scattering [ 22 ], or photopion 
production predominantly through the resonant process 
P 7 ^ A+ ^ n 7 r+ [17]. For the first of these, the in¬ 
elasticity of the process is 0.6 [23]. This then implies 
that the energy carried away by charged pions is about 
equal to the emerging nucleon energy, yielding (with our 
definition) e^r ~ 1. For resonant photoproduction, the in¬ 
elasticity is kinematically determined by requiring equal 
boosts for the decay products of the A+ [24], giving 

_I 


Ctt = E-K+lEn ~ 0.28, where if^+, E^ are the emerg¬ 
ing charged pion and neutron energies, respectively. 

In this section, we will extend the WB analysis [17] 
to the case of a lower onset of the extra-galactic com¬ 
ponent. The present analysis differs from WB in that 
the integrated power spectrum has changed, and that 
the spectral index 7 2. We will restrict the ensuing 

discussion to the case of photopion production on reso¬ 
nance, and comment on the pp possibility at the end of 
the section. 

The intermediate state of the reaction p -|-7 ^ fV-|- 7 r is 
dominated by the A"*" resonance. In order to normalize to 
the observed cosmic rays, we restrict our interest to the 
decay channel. Each tt’*' decays to 3 neutrinos and a 
positron, 7 r+ ^ . The e’*' readily loses 

its energy through synchrotron radiation in the source 
magnetic fields. The average neutrino energy from the 
direct pion decay is = (1 — r)if ^/2 ~ 0.22 

and that of the muon is {E^)t^ = (I -l-r) if^/2 ~ 0.78 
where r is the ratio of muon to the pion mass squared. 
Now, taking the from muon decay to have 1/3 the 
energy of the muon, the average energy of the from 
muon decay is (F^^)^ = (!-(- r)if^/6 = 0.26This 
means that neutrinos carry away about 3/4 of the 
energy, and each neutrino carries a fraction e^/4 of the 
accompanying cosmic ray energy. 

In order to correlate the neutrino and cosmic ray fluxes, 
we assume both follow a common power law at injection 


m 

dEi 


= C, 


CR(i^; 


Er 


( 4 ) 


and normalize our spectrum in a bolometric fashion 


^€7r£^2/4 q pE2 

c. / E-^-^-^Ue, = -e^CcR / . 

After integration we have. 






Gcr . 


( 5 ) 


Therefore, for the “low” crossover energy scenario, the 
resulting flux of neutrinos (all flavors) from optically thin 
sources is given by [17] 


ME) = 3 (I) 


^ ^ PCR dNn 


47r dEi 


Ei=E 


(I -^2)(«-7) 


dz ~ 3.5 X 10 


-3 


-^) 

GeV/ 


-2.54 


GeV ^cm ^s 


isr-b 


( 6 ) 


where the numerical value is an average flux obtained 
from best fits to AGASA and HiRes data, derived in 
Sec. II [25]. This extends the WB analysis to values of 
7 ^ 2 . The fluxes for each of the best fits are shown in 


Fig. 3, along with the WB flux. Also shown is the region 
excluded by AMANDA-BIO for both the cases 7 = 2 
and 7 = 2.54, and the cascade limit [26] from Ref. [27], 
which applies to all scenarios where neutrinos originate 
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from pion decays [28]. For the low-crossover scenario, the 
neutrino flux associated with the AGASA data set is con¬ 
sistent with the AMANDA-BIO data. Sensitivity to this 
flux will be attained by a full analysis of the AMANDA- 
II data set [30]. In the event of a positive indication 
by AMANDA-II, the IceCube facility will (as shown in 
Sec. V) be capable discriminating between the low and 
high crossover scenarios. 

Several other remarks are in order, (a) The best fit 
source neutrino fluxes in Fig. 3 correspond to Zmin = 
0.012 and GeV. We have checked that a 

change of Zmin to zero and a variation of within the 

range [10^^'^ GeV, 10^^'® GeV] produces changes within 
the order of the thickness of the lines. (5) Because of 
oscillations, the neutrino flavor mix at the source will 
evolve to a 1:1:1 mix at the detector [31]. (c) Electron 
antineutrinos can also be produced through neutron (3- 
decay. However, this contribution turns out to be neg¬ 
ligible (about 3 orders of magnitude smaller than the 
charged pion contribution) [32]. (d) The rapid rise at 

lower energies of the low crossover neutrino spectrum will 
greatly increase the event rate at the Glashow resonance 
as compared to that given in [32] based on the WB flux. 
The Vf. flux at the Glashow resonance obtained in the low 
crossover scenario is about an order of magnitude smaller 
than the bound from AMANDA-II data [33]. (e) Neu¬ 
trino fluxes from nearby isolated sources {e.g., Centau- 
rus A [22, 34] or Gygnus-OB2 [35]) can be distinguished 
through their point anisotropies, from the differential flux 
derived in this work. 

We now comment on the pp scenario. If the pri¬ 
mary proton spectrum oc E”"’’, the dominance of inelas¬ 
tic pp collisions produces an isotropically neutral mix of 
pions that on decay give rise to a neutrino flux with 
spectrum cx E~^ [36]. Gurrent hadronic event gener¬ 
ators yield an inelasticity of ~ 0.6 [23] for the reaction 
pp W-|-pions, where the V’s are final state nucleons. 
With our definition, « 1, assuming that 2/3 of the fi¬ 
nal state pions are charged. Then, the correction due to 
a larger inelasticity of pp interactions as compared to the 
resonant py scattering (with « 0.28) would increase 
the neutrino flux predictions given in Fig. 3 by a factor 
of « 7. This sizeable augmentation of the neutrino flux 
based on optically thin sources with dominant pp scatter¬ 
ing will result in the exclusion of the low crossover sce¬ 
nario. Therefore, we will continue to present our results 
on the basis of the dominance of the photopion process. 


IV. COSMOGENIC NEUTRINOS 

The opacity of the CMB to ultra-high energy pro¬ 
tons propagating over cosmological distances guarantees 
a cosmogenic flux of neutrinos, originated in the reac¬ 
tion p -I- 7CMB N -\- TT [37]. Very recently, one of us 
has performed an investigation of the actual size of the 
cosmogenic neutrino flux [13] assuming that all observed 
cosmic ray showers above 10®'^ GeV are initiated by pro- 
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FIG. 3: Neutrino fluxes (summed over all flavors) from 
optically thin sources for Ctt = 0.28. The horizontal solid 
line indicates the WB prediction which corresponds to a 
Galactic/extra-galactic crossover energy at the ankle, ~ 10^° 
GeV. The falling solid lines indicate the expected neutrino 
flux normalized to HiRes (lower) and AGASA (upper) data, 
if one assumes the onset of dominance by the extra-galactic 
component is at 10®’® GeV. The dash-dotted lines indicate 
the fluxes of cosmogenic neutrinos associated with flux pre¬ 
dictions given by the falling solid lines. The cross-hatched 
region excludes an E~^ spectrum at the 90% GL by measure¬ 
ments of AMANDA-BIO [29]. The single hatched region, ob¬ 
tained by rescaling the AMANDA integrated bolometric flux 
limit to an power law, is the exclusion region for the 

low crossover model. The shaded region indicates the cascade 
limit (see text for details). 


tons, with sources isotropically distributed throughout 
the universe. The low energy cutoff used in [13] is near 
the crossover energy suggested by HiRes data, and thus 
we expect no significant modification on the prediction 
of cosmogenic neutrinos. To verify this assertion, we es¬ 
timate the flux of neutrinos produced as sub-products 
in the GZK chain reaction by the population of protons 
that best reproduces the HiRes data. Such a cosmogenic 
flux is obtained by replacing Pp|„ in Eq. (1) with 


J. 



dP^\n{E-,Ei,r) 

dE 


( 7 ) 


In Eig. 4 it is seen that the cosmogenic flux predictions 
for the low energy crossovers at 10®'^ GeV and 10® ® GeV 
are compatible within errors. It should be noted that the 
contribution to the cosmogenic neutrino flux resulting 
from neutron beta decay is negligible for the energies 
under consideration [38]. 

The cosmogenic neutrino flux corresponding to the 
standard E~^ injection spectrum has been previously 
obtained in Ref. [38] for several assumed source evolu¬ 
tion indices. A comparison with the WB flux shows that 
these (the cosmogenic and source fluxes) are comparable 
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FIG. 4: Best fit to HiRes data (solid), assuming dominance 
of the extra-galactic component above E- = 10®'^ GeV (top) 
and E- — 10®'® GeV (bottom). Also sho'wn is the associated 
cosmogenic neutrino flux for all flavors (dashed). 


at energies above 10® GeV. This is in striking contrast 
■with the fluxes resulting from the low crossover scenario. 
As can be seen in Fig. 3, the source flux dominates the 
cosmogenic flux at energies below 10® GeV in the low 
crossover scenario [39]. Thus, the neutrinos below this 
energy behave as “unscathed messengers” of the source 
injection spectrum. The observation of a neutrino flux 
with a power law spectral index > 2.4 can provide strong 
support for the low crossover scenario. 


V. ICECUBE SENSITIVITY 

In the previous sections we have shown that if the 
nucleon-emitting sources are optically thin, then the dif¬ 
fuse flux of neutrinos produced by these sources provides 
a powerful tool in discriminating between Galactic/extra- 
galactic cosmic ray origin. Should the entire scenario not 
be ruled out by AMANDA-II data, it is of interest to ex¬ 
plore the potential of forthcoming neutrino telescopes to 
provide conclusive identification of the crossover energy. 

In deep ice/water/salt, neutrinos are detected by ob¬ 
servation of the Gerenkov light emitted by charged par¬ 
ticles produced in neutrino interactions. In the case of 
an incident high-energy muon neutrino, for instance, the 
neutrino interacts with a hydrogen or oxygen nucleus in 
the deep ocean water (or ice) and produces a muon trav¬ 
elling in nearly the same direction as the neutrino. The 
blue Cerenkov light emitted along the muon’s kilometer- 
long trajectory is detected by strings of photomultiplier 



E [GeV] 



E [GeV] 

FIG. 5: Upper panel: Differential event rate at IceGube for 
the different neutrino flnx predictions from AGASA (top), 
HiRes (middle) and WB (bottom) obtained in Sec. III. Lower 
panel: Expected bin-by-bin event rates for 10 years of opera¬ 
tion. The bin partition interval is taken as A logj^Q E — 0.5. 


tubes deployed at depth shielded from radiation. The 
orientation of the Cerenkov cone reveals the neutrino di¬ 
rection. There may also be a visible hadronic shower if 
the neutrino is of sufficient energy. 

The Antarctic Muon And Neutrino Detector Array 
(AMANDA) [40], using natural 1 mile deep Antarctic 
ice as a Cerenkov detector, has operated for more than 
3 years in its final configuration of 680 optical modules 
on 19 strings. The detector is in steady operation col¬ 
lecting roughly four neutrinos per day using fast on-line 
analysis software. Its performance has been calibrated 
by reconstructing muons produced by atmospheric muon 
neutrinos [41]. 

Overall, AMANDA represents a proof of concept for 
the kilometer-scale neutrino observatory, IceCube [18], 
now under construction. IceCube will consist of 80 
kilometer-length strings, each instrumented with 60 10- 
inch photomultipliers spaced by 17 m. The deepest mod¬ 
ule is 2.4 km below the surface. The strings are arranged 
at the apexes of equilateral triangles 125 m on a side. The 
instrumented (not effective!) detector volume is a cubic 
kilometer. A surface air shower detector, IceTop, consist¬ 
ing of 160 Auger-style [42] Cerenkov detectors deployed 
over 1 km® above IceCube, augments the deep-ice com- 
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ponent by providing a tool for calibration, background 
rejection and air-shower physics. Construction of the de¬ 
tector started in the Austral summer of 2004/2005 and 
will continue for 6 years, possibly less. At the time of 
writing, data collection by the first string has begun. 

At the energies under consideration, there is no atmo¬ 
spheric muon or neutrino background in a km^ detector. 
The differential event rate is given by 

~ 27rACi p Veff Ji, ! (8) 

where Na is Avogadro’s number, Veff « 2 km^ is the ef¬ 
fective volume of ice with density p, and = 6.78 x 
10“^® (i<l/TeV)°'^®^ cm^ is the charged current neutrino- 
nucleon cross section [43]. The effective volume used is 
conservative, since muon tracks can originate well out¬ 
side the fiducial volume of the detector [44]. In Fig. 5 
we show the differential event rate at IceCube from op¬ 
tically thin sources. Also shown are the expected bin- 
by-bin event rates for 10 years of data collection, with 
a bin partition size Alogj^pA = 0.5. The vertical error 
bars are obtained on the basis of Poisson statistics with 
A logj^Q N = 0.434-\/]V /N, for N > 20. For smaller statis¬ 
tics we use Poisson confidence intervals [50]. It is strongly 
indicated that within its lifetime IceCube will attain suf¬ 
ficient sensitivity to constrain the energy of transition be¬ 
tween Galactic and extra-galactic dominance. RICE [45], 
PAO [46], EUSO [47], ANITA [48], and OWL [49] also 
have the potential to measure the ultra-high energy neu¬ 
trino flux. However, the energy thresholds, systematics, 
backgrounds, or time-scales to completion leave these ex¬ 
periments less promising than IceCube for a spectrum 
determination in the near future. 


VI. CONCLUSIONS 

We have estimated the extragalactic diffuse neutrino 
flux emitted from optically thin sources, on the basis 
of a low transition energy (10® ® GeV) between Galactic 
and extragalactic cosmic rays. Such a low crossover finds 
support in the chemical composition analysis of HiRes 
data [7], and is sustained by studies which reproduce the 
steepening at the second knee via e+e“ production on 
the CMB [11]. Since the neutrino flux reflects the nucleon 
flux at the source, the latter must be obtained by fitting 


observed cosmic ray data taking into account propaga¬ 
tion effects. The low crossover energy is well below the 
threshold energy for resonant pjcMB absorption, and so 
samples sources even at large redshift. Thus, source evo¬ 
lution is an important consideration in this calculation. 
Two-parameter fits in the spectral and redshift evolu¬ 
tion indices show that a standard Eermi E~'^ source spec¬ 
trum is excluded at larger than 95% CL. Best fits to both 
Akeno [9] -I- AGASA [2] and Fly’s Eye [8, 15] -I- HiRes [4] 
data sets give an source spectrum, with an evo¬ 

lution index somewhat larger than 3. The neutrino flux 
obtained using the WB [17] consideration for energetics 
at the source mirrors the steep spectrum of the emitted 
cosmic rays. 

Comparison of the resulting flux with existing 
AMANDA-BIO 90% CL bounds [29] reveals the follow¬ 
ing: (1) If neutrinos are generated by pp interactions at 
the source, the resulting flux is within the excluded re¬ 
gion. (2) For pj interactions dominant, the best fit to 
the data yields a neutrino flux which is consistent with 
the AMANDA-BIO upper limit. A complete analysis of 
the AMANDA data will provide sufficient sensitivity to 
rule out the model. 

The neutrino flux at the source in this scenario dom¬ 
inates the cosmogenic flux. Thus, should data from 
AMANDA-H not rule out the model, we show that Ice- 
Cube can measure the characteristic power law of the 
neutrino spectrum, and thus provide a window on the 
source dynamics. 

In summary, forthcoming data from the South Pole 
can provide significant clues in demarcating the cosmic 
ray Galactic extra-galactic crossover energy. 


Acknowledgments 

We would like to thank Birgit Eberle and Yvonne 
Wong for discussions, and Zoltan Fodor and Sandor Katz 
for supporting our numerical analysis based on their orig¬ 
inal computer codes. This work has been partially sup¬ 
ported by US NSF (Grant Nos. OPP-0236449, PHY- 
0140407, PHY-0244507), US DoE (Grant Nos. DE- 
EG02-95ER40896, DE-EG05-85ER40226), NASA-ATP 
02-000-0151, and the Wisconsin Alumni Research Eoun- 
dation. 


[1] See e.g., L. Anchordoqui, T. Paul, S. Reucroft and 
J. Swain, Int. J. Mod. Phys. A 18, 2229 (2003) 
[arXiv:hep-ph/0206072]. 

[2] M. Takeda et al, Phys. Rev. Lett. 81, 1163 

(1998) [arXiv:astro-ph/9807193]; M. Takeda et al., As- 
tropart. Phys. 19, 447 (2003) [arXiv:astro-ph/0209422]. 
AGASA (Akeno Giant Air Shower Array) http://www- 
akeno.icrr.u-tokyo.ac.jp/AGASA/ 24 February 2003 


[3] K. Greisen, Phys. Rev. Lett. 16, 748 (1966); G. T. Zat¬ 
sepin and V. A. Kuzmin, JETP Lett. 4, 78 (1966) [Pisma 
Zh. Eksp. Teor. Fiz. 4, 114 (1966)]. 

[4] T. Abu-Zayyad et al. [HiRes Gollaboration], arXiv:astro- 
ph/0208301. 

[5] N. Hayashida et al. [AGASA Collaboration], Phys. Rev. 
Lett. 77, 1000 (1996); N. Hayashida et. al. [AGASA 
Collaboration], arXiv:astro-ph/0008102; G. R. Farrar, 



arXiv;astro-ph/0501388. See, however, C. B. Finley 
and S. Westerhoff, Astropart. Phys. 21, 359 (2004) 
[arXiv: astro-ph/0309159]. 

[6] R. U. Abbasi et al. [HiRes Collaboration], arXiviastro- 
ph/0404137; R. U. Abbasi et al. [HiRes Collaboration], 
arXiv:astro-ph/0404366. R. U. Abbasi et al. [HiRes Col¬ 
laboration], arXiv:astro-ph/0412617. 

[7] D. R. Bergman [HiRes Collaboration], astro-ph/0407244. 

[8] D. J. Bird et al. [Fly’s Eye Collaboration], Phys. Rev. 
Lett. 71 3401 (1993). 

[9] M. Nagano et al., J. Phys. G 18, 423 (1992). 

[10] T. Abu-Zayyad et al., [HiRes-MIA Collaboration], Astro- 
phys. J. 557, 686 (2001) [arXiv:astro-ph/0010652]. 

[11] V. Berezinsky, A. Z. Gazizov and S. I. Grigorieva, 
arXiv:hep-ph /0204357. 

[12] R. U. Abbasi et al. [HiRes Gollaboration], Phys. 

Rev. Lett. 92, 151101 (2004) [arXiv:astro-ph/0208243[; 
V. Berezinsky, A. Gazizov and S. Grigorieva, astro- 
ph/0302483; V. S. Berezinsky, S. 1. Grigorieva and 
B. 1. Hnatyk, Astropart. Phys. 21, 617 (2004) 

[arXiv:astro-ph/0403477]; V. Berezinsky, A. Z. Gazizov 
and S. 1. Grigorieva, arXiv:astro-ph/0502550. 

[13] Z. Fodor, S. D. Katz, A. Ringwald and H. Tu, JCAP 
0311, 015 (2003) [arXiv:hep-ph/0309171]. 

[14] D. Seckel and T. Stanev, arXiv:astro-ph/0502244. 

[ 15 ] D. J. Bird et al. [HiRes Gollaboration], Astrophys. J. 424, 
491 (1994); D. J. Bird et al., Astrophys. J. 441, 144 
(1995). 

[16] Z. Fodor and S. D. Katz, Phys. Rev. D 63, 023002 
(2001) [arXiv:hep-ph/0007158]; Z. Fodor, S. D. Katz 
and A. Ringwald, JHEP 0206, 046 (2002) [arXiv:hep- 
ph/0203198]; Z. Fodor, S. D. Katz, A. Ringwald and 
H. Tu, Phys. Lett. B 561, 191 (2003) [arXiv:hep- 
ph/0303080]; http://www.desy.de/'uhecr 

[17] E. Waxman and J. N. Bahcall, Phys. Rev. D 59, 023002 
(1999) [arXiv:hep-ph/9807282]. 

[18] J. Ahrens et al. [IceGube Collaboration], Astropart. 
Phys. 20, 507 (2004) [arXiv:astro-ph/0305196]. 

[19] E. Waxman, Astrophys. J. 452, LI (1995). See also, 
T. K. Gaisser, AIP Conf. Proc. 566, 238 (2000) 
[arXiv:astro-ph/OOl 1525]. 

[20] A. Miicke, R. Engel, J. P. Rachen, R. J. Protheroe and 
T. Stanev, Comput. Phys. Commun. 124, 290 (2000) 
[arXiv:astro-ph/9903478]. 

[21] D. N. Spergel et al. [WMAP Collaboration], Astrophys. 
J. Suppl. 148, 175 (2003) [arXiv:astro-ph/0302209]. 

[22] L. A. Anchordoqui, H. Goldberg, F. Halzen and 
T. J. Weiler, Phys. Lett. B 600, 202 (2004). 

[23] J. Alvarez-Muniz, R. Engel, T. K. Gaisser, J. A. Ortiz 
and T. Stanev, Phys. Rev. D 66, 033011 (2002). 

[24] F. W. Stecker, Phys. Rev. Lett. 21, 1016 (1968). 

[25] The derivation of the kernel of integration in Eq. (6) can 
be obtained using the procedure presented by B. Eberle, 
A. Ringwald, L. Song and T. J. Weiler, Phys. Rev. D 70, 
023007 (2004) [arXiv:hep-ph/0401203]. 

[26] V. S. Berezinsky and A. Yu. Smirnov, Astrophysics and 
Space Science 32, 461 (1975). 

[27] K. Mannheim, R. J. Protheroe and J. P. Rachen, Phys. 
Rev. D 63, 023003 (2001) [arXiv:astro-ph/9812398]. 

[28] The cascade limit arises from the requirement that the 
diffuse gamma-ray fluxes initiated by photons and elec¬ 
trons from pion decays should not exceed measure¬ 
ments. The cascade limit from Ref. [27] exploits the mea¬ 


surement of the diffuse 7 ray background by EGRET, 
P. Sreekumar et al., Astrophys. J. 494 (1998) 523 
[arXiv:astro-ph/9709257]. A lower extragalactic contri¬ 
bution to the 7 ray background than that inferred in 
abover reference, by roughly a factor of two, has been 
proposed recently by A. W. Strong, 1. V. Moskalenko and 
O. Reimer, arXiv:astro-ph/0306345; U. Keshet, E. Wax- 
man and A. Loeb, JGAP 0404, 006 (2004) [arXiv:astro- 
ph/0306442]. The cascade limit may therefore be stronger 
by a corresponding factor. 

[29] M. Ackermann et al., Astropart. Phys. 22, 339 (2005). 

[30] F. Halzen, talk presented at XI International Workshop 
on Neutrino Telescopes, February 2005. 

[31] For alternative neutrino oscillation models see e.g., 
J. F. Beacom, N. F. Bell, D. Hooper, S. Pakvasa and 
T. J. Weiler, Phys. Rev. D 68 , 093005 (2003) [arXivihep- 
ph/0307025]. 

[32] L. A. Anchordoqui, H. Goldberg, F. Halzen and 

T. J. Weiler, arXiv:hep-ph/0410003. 

[33] M. Ackermann [AMANDA Collaboration], Astropart. 
Phys. 22 , 127 (2004) [arXiv:astro-ph/0405218]. 

[34] L. A. Anchordoqui, H. Goldberg and T. J. Weiler, Phys. 
Rev. Lett. 87, 081101 (2001) [arXiv:astro-ph/0103043]. 

[35] L. A. Anchordoqui, H. Goldberg, F. Halzen and 

T. J. Weiler, Phys. Lett. B 593, 42 (2004) [arXiv:astro- 
ph/0311002]. 

[36] We stress that logarithmic corrections to the spectrum 
oc E ^ could result in sizeable deviations at very high 
energy. L. Anchordoqui, H. Goldberg and C. Nunez, 
arXiv:hep-ph/0408284. 

[37] V. S. Berezinsky and G. T. Zatsepin, Phys. Lett. B 28 
423 (1969). 

[38] R. Engel, D. Seckel and T. Stanev, Phys. Rev. D 64, 
093010 (2001) [arXiv:astro-ph/0101216]. 

[39] We have checked that for a change of Zmin to zero and 

for larger than the default 10^^ ® GeV, the result¬ 

ing cosmogenic neutrino fluxes do not change within the 
energy range shown in Fig. 3. For smaller Ei^max, the 
predicted fluxes get even smaller. 

[40] E. Andres et al. [AMANDA Collaboration], Astropart. 
Phys. 13, 1 (2000) [arXiv:astro-ph/9906203]. 

[41] E. Andres et al., Nature 410, 441 (2001). 

[42] J. Abraham et al. [Auger Collaboration], Nucl. Instrum. 
Meth. A 523, 50 (2004). 

[43] R. Gandhi, C. Quigg, M. H. Reno and 1. Sarcevic, Phys. 
Rev. D 58, 093009 (1998) [arXiv:hep-ph/9807264]. 

[44] F. Halzen and D. Hooper, Rept. Prog. Phys. 65, 1025 
(2002) [arXiv:astro-ph/0204527]. 

[45] 1. Kravchenko et al. [RICE Collaboration], Astropart. 
Phys. 19, 15 (2003) [arXiv:astro-ph/0112372]. 

[46] K. S. Capelle, J. W. Cronin, G. Parente and E. Zas, 
Astropart. Phys. 8 , 321 (1998) [arXiv:astro-ph/9801313]. 

[47] L. Scarsi, Nuovo Cim. 24C, 471 (2001); O. Catalano, 
Nuovo Cim. 24C, 445 (2001). 

[48] A. Silvestri et al. [ANITA Collaboration], arXiv:astro- 
ph/0411007. 

[49] F. W. Stecker, J. F. Krizmanic, L. M. Barbier, E. Loh, 
J. W. Mitchell, P. Sokolsky and R. E. Streitmatter, 
Nucl. Phys. Proc. Suppl. 136C, 433 (2004) [arXiv:astro- 
ph/0408162]. 

[50] G. J. Feldman and R. D. Cousins, Phys. Rev. D 57, 3873 
(1998) [arXiv:physics/9711021]. 



